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CALCULATED PRESSURE-INDUCED VIBRATIONAL ABSORPTION 
IN h 2 - h 2 COLLISIONS IN HYDROGEN GAS 
by R. W. Patch 
Lewis Research Center 

SUMMARY 

The coefficient for pressure- induced vibrational absorption in H 2 - H 2 collisions 
was calculated for temperatures from 298 to 7000 K and wave numbers between 100 and 
40 000 reciprocal centimeters for local thermodynamic equilibrium. Because only 
transitions with a net change of +1 vibrational quanta were considered, the absorption 
was centered near the fundamental at 4161 reciprocal centimeters in the infrared. The 
model included electronic configuration interaction, mechanical anharmonicity, 
vibration- rotation interaction, excited vibrational states, and more realistic intermo- 
lecular potential and line shapes than previously used. The integrated absorption coeffi- 
cient at 298 K was 87 percent of the experimental value. The integrated absorption coef- 
ficient at 3000 K was 2. 1 times the previous theoretical value. For pressures from 100 
to 1000 atmospheres (1. 013x10^ to 1. 013x10^ N/m^) the Planck and Rosseland mean 
opacities were increased by factors as high as 18. 9 and 303, respectively, over pre- 
vious estimates. An approximate formula for the spectral absorption coefficient is given 
for rapid calculation for temperatures from 600 to 7000 K. 


INTRODUCTION 

The need for reliable heat-transfer calculations for high- temperature high-pressure 
hydrogen gas occurs in gaseous-core nuclear rockets (refs. 1 and 2), in outer layers 
of late-type stars, and in high-speed entry of spacecraft into certain planetary atmo- 
spheres. Such calculations require accurate pressure-induced vibrational absorption 
coefficients for H 2 - H 2 collisions. This absorption occurs in the infrared region and is 
principally due to transitions with a net change of +1 vibrational quanta for the two mole- 
cules. The special case where the vibrational quantum number of one molecule jumps 
from 0 to 1 is known as fundamental vibrational absorption. 
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Considerable theoretical work has been performed on pressure -induced Hg - Hg 
fundamental vibrational absorption. Isolated Hg molecules do not possess an electric 
dipole moment, but a necessary condition for pressure-induced vibrational absorption is 
an electric dipole moment. This moment is induced during collisions. Van Kranendonk 
and Bird (ref. 3) gave a theory which assumed that the induced dipole moment was the 
sum of a contribution due to the overlap of orbitals and a contribution due to the electric 
quadrupole moment of each molecule inducing a dipole moment in the other molecule. 

This is called the large-R approximation in this report. They evaluated the overlap con- 
tribution by a distorted- orbital method. Britton and Crawford (ref. 4) improved Van 
Kranendonk and Bird's theory by including components of the dipole moment perpendicu- 
lar to the intermolecular axis and by including an additional exchange effect in calculating 
the overlap dipole moment. Van Kranendonk (refs. 5 and 6) extended his theory to cryo- 
genic temperatures by considering quantum effects in the translational motion of the gas 
molecules. Linsky (ref. 7) extrapolated Van Kranendonk's (ref. 6) theory to tempera- 
tures up to 3000 K although the theory neglected excited vibrational states. All theories 
have assumed a Lennard-Jones 6-12 intermolecular potential (ref. 8). 

A great deal of experimental work has also been performed on pressure -induced 
Hg - fundamental vibrational absorption. Experimental measurements of absorption 
coefficients have been made at temperatures up to 423 K (refs. 9 to 11) and are in gen- 
eral agreement with each other and with theory (refs. 4 and 6). 

At temperatures several times higher than existing experimental data, several ef- 
fects can be expected to be important that were not important in the experiments or 
theories. Excited vibrational states may be more important than their relative popula- 
tions would indicate, and mechanical anharmonicity and vibration- rotation interaction 
should also have appreciable effects. In addition, as temperature is raised, collisions 
become harder so that induced dipole moment and intermolecular potential must be 
known for smaller intermolecular distances. It is not safe to extrapolate dipole moment 
from values for intermolecular distances important at room temperature because elec- 
tronic configuration interaction is important at small intermolecular distances but not at 
large intermolecular distances. It is well known not to be safe to use a Lennard-Jones 
6-12 intermolecular potential at high temperature because it gives excessive potentials 
for small intermolecular distances. Lastly, at high temperature the absorption coeffi- 
cient must be reduced because of the effect of stimulated emission, which follows 
Kirchhoff's law. This has not been done correctly in some previous theories. 

The object of this report was to calculate more realistic values of absorption coeffi- 
cient for transitions with a net change of +1 vibrational quanta for temperatures up to 
7000 K. All the above high-temperature effects were included as realistically as cur- 
rently possible. The induced dipole moment was obtained from orthogonaliz ed valence- 
bond configuration-interaction calculations for small intermolecular distances (ref. 12) 
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and by the large-R approximation for large inter molecular distances (ref. 13). The ef- 
fects of excited vibrational states, mechanical anharmonicity, and vibration-rotation in- 
teraction were based on a new analysis presented in this report. (Configuration inter- 
action, excited vibrational states, mechanically anharmonicity, and vibration- rotation 
interaction were not included in previous pressure-induced absorption theories. ) More 
recent line shapes were used than Linsky (ref. 7) used. The temperature range 298 to 
7000 K was covered. At higher temperatures other processes are much more important. 
This report is limited to densities where ternary collisions are not important (ref. 10). 


ANALYSIS 

In this section the integrated absorption coefficient is obtained for pressure- induced 
vibrational transitions including the effects of excited vibrational states, mechanical an- 
harmonicity, and vibration- rotation interaction. The relation of this coefficient to sev- 
eral other absorption coefficients is then given, and realistic line shapes are selected 
based on experiments in the literature and on theoretical considerations. 


Absorption Coefficients in Terms of Einstein Coefficients 


The first step in getting the integrated absorption coefficient in useful form is to re- 
late the linear absorption coefficient, the Einstein coefficients for absorption, and the in- 
tegrated absorption coefficient. Consider a gas containing only Hg molecules. We are in- 
terested in the pressure-induced vibrational absorption occurring during binary collisions 
of pairs of these molecules. Let v^, J^, and rm be the vibrational quantum number, total 
angular momentum quantum number, and quantum number for the component of total an- 
gular momentum along the inter molecular axis, respectively, for the i^ molecule. (Sym- 
bols are given in the appendix. ) A state p of the pair thus has the quantum numbers v^, 
Vg, Jj, Jg, m p and mg provided the molecules are far enough apart. We define a spec- 
tral Einstein coefficient of absorption B in such a way that the probability of a pair of 
molecules in state p, exposed to radiation of wave number v, absorbing a quantum hep 
in time interval dt and wave number interval dv and making a transition to state p' is 
given by B^,~p~ dt dv, where p~ is the radiation energy density per unit wave number. 

According to Chandrasekhar (ref. 14), the linear absorption coefficient a~ (excluding 
stimulated emission) is related to the spectral Einstein coefficient by 


a~ = p B — 
v *p pp'v y 


( 1 ) 


3 



where V is the volume of homogeneous gas considered, is the number of pairs in 
state p, and h is Planck's constant. The integrated Einstein coefficient of absorption 
for the transition p — p' is 



( 2 ) 


The integrated absorption coefficient with the wave number removed is defined by 


pp' 



Combining equations (1) to (3), 


P hB , 
s - P PP 

V 


’PP’ 


(3) 


(4) 


However, B , actually depends on the inter molecular distance R of the pair of mole- 
cules (fig. 1), so equation (4) must be corrected by averaging over R. For a dilute sys- 
tem the averaging may be accomplished by (ref. 15) 


P h 

s = _P 
PP' 2 
V z 


/ 


B e “ $//kT 47rR 2 dR 
PP 


(5) 


where $ is the intermolecular potential averaged over orientations so that it depends 
only on R. For neutral molecules $ decreases rapidly with increasing R, so the inte- 
gration may be extended to R = 00 without loss of accuracy. 

More useful versions of equation (5) are obtained by replacing by number densi- 
ties and summing. The number of pairs of H 9 molecules is N(N - l)/2, where N is the 

^ 2 

number of such molecules in volume V. For N » 1 this is essentially N /2. The 
quantity is thus 


r N 


p p = lT/ P(Vl ’ Jl ’ mi)P(V2,J2, m2) 


( 6 ) 
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where P(v, J, m) is the probability of an Hg molecule having quantum numbers with val- 
ues v, J, and m. States with the same Vp Vg, Jp and Jg but any values of m^ 
and mg all have the same energy. Likewise, states with the same vj, Vg, Jp and Jg 
but any values of m j and mg all have the same energy. Consequently, transitions be- 
tween these primed and unprimed states all have the same wave number. Thus for con- 
venience their integrated absorption coefficients (eq. (5)) may be added. 


s zv 




J 2 J 2 

Y Y P(V 1’ J V m l> P(v 2’ J 2’ m 2> 

m 2 =-J 2 m 2 =-J 2 


< 7 > 

where subscript £ stands for the quantum numbers Vp Jp Vg, and Jg; a prime indi- 
cates quantum numbers after the transition: and n is the number density corresponding 
to N. 


Wave Functions and Expansions 

The Einstein coefficient and hence s ^> may be expressed in terms of the wave 
functions and of expansion coefficients of the dipole moment. First a set of coordinates 
are needed. We neglect translation of the molecules and choose Cartesian coordinates 
fixed in space half way between the molecules (fig. 1). 

The integrated Einstein coefficient is related to the components ( n \ , (ii \ , 

z x ' X/ P'P ' y 'P'p 

and ( p ) °f the electric dipole matrix elements for these coordinates by (ref. 16) 

V z Vp 



(8) 


where e Q is the electric permittivity of free space, c is the velocity of light, and the 
equation is in SI units. The x and y contributions to equation (8) can be expected to be 
considerably smaller than the z contribution for two reasons: (1) The electric quadru- 
ple moment of each molecule may induce a dipole in the other molecule. The average 
magnitudes of the x and y components of the electric fields due to these quadrupole 
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moments are much less than the z components. (2) The distortion of the orbitals due 
to overlap is principally in the z direction. Therefore the x and y contributions to 
equation (8) are neglected here. If equation (8) is substituted into equation (7), we get 


s rr 


A 2 P( Vl ,J r 0)P(v 2 ,J 2 , 0) 

3hC€^ 

o 





x /n\ 2 e ~ 9//kT 4ttR 2 dR (9) 

V Jp'p 


where the fact that P(v,J, m) is independent of m has been used. 

The relation between the matrix element ( p ) and the system wave functions \p 

a , ■ ' ' p’p P 

and \l/ , is H K 

V 



th th 

where z. is the z coordinate of the l electron, z- is the z coordinate of the j m 

1 J th 

nucleus, e is the charge of an electron, eZj is the charge of the j nucleus, and the 

summations cover all electrons and nuclei. If the relative motion of the two molecules 

is slow compared to the motion of the electrons, the Born-Oppenheimer approximation 

(ref. 17) may be applied to the pair of molecules collectively. 


= * e (X., R, r 1? r 2 , <p v e 2 ,(p 2 )\l/ n (R,r 1 ,r 2 , <p v 6> 2 , cp 2 ) 


( 11 ) 


where \Is is the electronic wave function, and x V is the nuclear wave function. In 
6 n 

equation (11) stands for the position coordinates of the electrons, r, for the inter- 
nuclear distance of molecule 1, and r 2 for the internuclear distance of molecule 2. It 
has been observed that the wave numbers of lines in pressure-induced transitions can be 
calculated from term values for isolated diatomic molecules. Consequently, it is rea- 
sonable to approximate \)/ n by (ref. 3) 


<p, 


Vi 1 *' 1 W 2 ’ 


n 


^l) Y J 2 m 2 ^2 ,<p 2^ 


( 12 ) 
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where i// y j is the vibrational wave function. In the past harmonic oscillator wave func- 
tions have been used for i// y j, but in this report mechanical anharmonicity is included 
by using a wave function calculated from the Rydberg- Klein- Rees potential energy for an 
isolated Hg molecule. In addition, vibration- rotation interaction in H2 is included by 
2 2 2 

adding h J(.J + 1)/8 t r m r r to the potential energy (ref. 18) before calculating the vibra- 
tional wave function. In equation (12) Yj m is the spherical harmonic given by 


Y Jm (e > V) = (2,)- 1/2 © T ™(cos e)e im< ? 


J nr 


(13) 


For m > 0 the normalized associated Legendre functions © Jm are given by Pauling 
and Wilson (ref. 19). For m < 0 they are 


©jm( c °s B) = (-l) m 0j I m l (cos 9) 


(14) 


From equations (10) to (12) 




P P 





X VlV2 Yj l m l Tj 2"2 Sln 61 S ‘ n 02 dri dr2 dSl d<, 2 d?, 2 <15) 


where the z component of the electric dipole moment is given by 


ju = 


J Z 6Z i + Z eZ i Z i 


^e dT e 


(16) 


Because z. is a linear, Hermitian operator, p must be real. The dipole moment a 
l z ~z 

may be expanded in a Taylor series about the equilibrium internuclear distances rV 

A 


and r 0 of the molecules. 


P-2^^’ r l’ r 2’ ^1’ *^l’ ^2’ M R > ^1’ ^ ®2’ ^2) + ^2’ < ? ? 2^ r 1 - l^ 


+ M 2 < R ’ e V Vv e 2’ ^2^ r 2 ■ r 2^ 


(17) 
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where p° is the value of u z for r° and r 2 , is 9M z / 9r i> and M 2 is 9 M z / 9r 2> 

both evaluated for r° and r^. The quantities p° and pj can be expanded in spherical 
harmonics. 


*1 _ '2 


m z (R, o v < p v 0 2 » < p $ ~ ^ ^ 2ffC *i*i l 2*s 


(R) 


l \~ 0 ^ 2 _0 ^ 2~~ l 2 


x Y 




ii l 2 

OO i- OO Z 


M 


j(R, «j, e 2 , v 2 ) = ^ ^ ^ 2 " r V 1 


^2^2 


(R) 


ij-0 ^ ^2 *2 


X 


Y l, ^ Y l 2 l 2 < 9 2' f 2 ) (19) 


where C 7 t 7 t and D 7 t 7 t are expansion coefficients to be determined. 

*in*2*2 t i«r2*2 


Molecular Symmetry and Configurations 

Considerable simplification of equations (18) to (19) is possible. Some of the expan- 
sion coefficients are zero because of molecular symmetry. The operation 


<?. - <?. + 7T 




e i 


which is equivalent to interchanging the two identical nuclei, should not change u or 


Mi 

cp 2 


for i = 1 or 2. In addition, p and p. should have azimuthal dependence solely on 

Z A 

- <p* because the orientation of the x and y axes should have no effect. Thus, the 


only nonzero C 7 t , t are 
4 i*r2«2 


' 222 - 2 ’ 


and Ci t 1 t with i-i — 4 or 

l ih l 2*2 1 


C 0000’ C 2000’ C 0020’ C 2020’ C 2-121» 

= C. 


^2 ~ 4 " 


C 212- 1’ C 2-222’ 
o 


Also, C 2 _i 2 i = C 2 i 2 _i because p 
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must not depend on what direction about the z axis cp ^ and <p 2 
larly C 2 _ 2 22 ~ ^222-2" ^11 the preceding findings apply to the 


are measured. Simi- 

D, f , , as well. 
Z 1«1 Z 2«2 


Further simplification of equations (18) to (19) is necessary because calculations of 
ju are time consuming even on a high-speed digital computer. This is especially true 
if small values of R are involved, which occurs at high temperature. Thus it is neces- 
sary to limit the number of configurations 9^, <7, # 2 , <?2 to as ^ ew as P oss ^t»le. The 
ones selected are shown in figure 2. The and ©2_j are zero for all these config- 

urations, so C 2-121’ ^212-1’ ^2-121’ anc * ^212 1 canno ^ be determined. There are 
just enough configurations to determine all remaining nonzero C 7 t , t and D, £ 7 t 

with < 4 and *= 4. Neglect of expansion coefficients with l j > 4 and/or ^ 4 
is reasonable because, if such coefficients were nonzero they would allow transitions 
with Jj or Jg changing by ±4 or more, which have never been observed. For r j = r° 
and rg = Tg, p z is zero for configurations 1, 2, and 5 by symmetry and differs only in 
sign for configurations 3 and 4. Hence equation (18) gives 


a 1 n , V5 P , Vs p , 5 p qp 

U “ 2 C ° 000 ~2 ° 2000 ~2 ° 0020 2 02020 2 " 222 

0 = - C — 


a p » ^ p Vs p , 5 

2 '-'0000 " y- zOOO " 


V5 C , 5 C 15 

7^0020 g 2020 g '"'2-222 


,.0 _ 1 p , Vs p V 5 p _ S p QP 

"^z, 4 “ 2 0000 7" 2000 4 ^0020 4 L 2020 2-222 

ii° - ^ p V5 p Vs p _ S p + qp 

^z, 4 “ “ C 0000 “ 7- u 2000 + ~ '-'0020 4 ^2020 + u '"2-222 

0 - 1 p Vs p V 5 p + S p 15 p 

U C 0000 “ ~ 7 ~ '"'2000 ~T ^0020 + A '-'2020 '-' 2-222 


( 20 ) 


This has the solution C 0000 = C 2020 = C 2 _ 222 = 0, C 2Q00 = -C 0Q20 . 

The derivative p 2 can be expressed in terms of pj. To do this, an arbitrary con- 
figuration is reflected in a plane bisecting the intermolecular axis, with the result 


p. 2 (H> <p 9 2 , < 5 ^ 2 ) ^ — ^2’ ^2’ ^ ” ^1’ 


( 21 ) 


where the order of the arguments of p^ has the significance given in equation (19). 
Substituting equation (21) into equation (17) gives 
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M z = M°( R > 6 i> <Pi> d 2> + Mj( R , 9 1> e 2’ ^2^ r l “ r °) 


- p^R, * - e 2 , <p 2 , n - 9 V ^j)(r 2 - r°) (22) 


Evaluation of Integrals and Sums 

The quadruple sum of the squares of the dipole matrix elements in equation (9) can 
be expressed simply in terms of vibrational overlap integrals, vibrational matrix ele- 
ments, expansion coefficients, and Kronecker delta functions. 

The first step is to combine equations (15), (18), (19), and (22) with the result 

K) = M 1 A 00^ A 00^ + M 2 A 20^ A 00^ + M 3 A 00^ A 20^ + M 4 A 20^ A 20^ 

\ /p’p 

+ ^5 A 2- 2^ ^ A 22^ + ^5 A 22^ ^ A 2- 2^ 2 ^ (23) 


where 


A ij( k) 


/ 7 T 


6 K’ j+m k )/ ^m^ij 0 ^ 


kk 


sin 0^ d0^ 


(24) 


M 1 3 ( A 1^2 " ^1 A 2^ D 0000 
M 2 3 < ^ 1 < ^ 2 C 2000 + A 1 ^ 2 D 2000 ' ^1 A 2 D 0020 
M 3 3 ~ Q 1 Q 2 C 2000 + A 1 Q 2 D 0020 ' Q 1 A 2 D 2000 
M 4 3 ^ A 1^2 “ ^1 A 2^ D 2020 

Mg s (AjQ 2 - Qj A 2 )D 2 222 (25) 


% ' r 

Jo 


^v ’ J 

Vk 


,\Lf T dr. 
% J k k 


( 26 ) 
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(27) 



^k J k (l k " r ^ 


Vk 


dr k 


The quantities and are the vibrational overlap integral and vibrational matrix 
element, respectively. The quantity 6(m£, j + m^) is the Kronecker delta function. The 
A„(k) are readily evaluated by using the properties (ref. 20) and orthogonality of the 
normalized associated Legendre functions. 


A 0 0 (k) = 2 1//2 6(m’, m)6(J',J) 


(28) 


A 20^ " 


/ r i 1/2 

Vio f KJ + l) 2 - m 2 l[(J + 2) 2 - m 2 l p T u R 2 

^ 3<^ 5(J’,J + 2) + M + 1)- 6m 5(J>,J) 

4 1 ] (2J + 1)(2J + 3) 2 (2J + 5) J + - X ) 


., f[ j 2 -m 2 ][(J-l) 2 -m 2 ] 
(2J + 1)(2J - 1) 2 (2J - 3) 


1/2 


6(J’,J - 2) 6(m', m) (29) 


A 22^ ~ 


_ Vl5 


(J + m + 4)(J + m + 3)(J + m + 2)(J + m + 1) 
(2J + 5)(2J + 3) 2 (2J + 1) 


1/2 


6(J',J + 2) 


- 2 


(J + m + 2)(J + m + 1)(J - m)(J - m - 1) 
(2J + 3) 2 (2J - l) 2 


1/2 

6(J’,J) 


(J - m)(J - m - 1)(J - m - 2)(J - m - 3) 


(2J + 1)(2J - 1)^(2J - 3) 


1/2 


6(J’, j - 2) 


>6(m’,m + 2) (30) 


1 1 



2-2 


(k) = 


Vl5 


(J - m + 4)(J - m + 3)(J - m + 2)(J - m + 1) 


(2J + 5)(2J + 3)^(2J + 1) 


1/2 


6(J', J + 2) 


- 2 


(J - m + 2)(J - m + 1)(J + m )(J + m - 1) 
(2J + 3) 2 (2J - l) 2 


11/2 


6(J’,J) 


(J + m)(J + m - 1)(J + m - 2)(J + m - 3) 


- 11/2 


(2J + 1)(2J - 1)^(2J - 3) 


6(J’,J- 2) U(m’, m- 2) 


(31) 


where the subscript k is understood to apply to all quantities on the right sides of equa- 
tions (28) to (31). 

Squaring and summing equation (23) and making use of equations (28) to (31) gives 


j ! j; j 9 Jo 

in ii <->„ 

m r- J l m r“ J i m 2 = - J 2 m 2=- J 2 


j, J i) L o< J 2t J 2> 


+ ^2^2^ 1 > ^ 1 ^ 0 ^ 2 ? ^ 2^ + M 3 V J i J ^ 1^2^ 2 7 ^ 2 ^ 


(m| + 2M^L 2 (J|,J 1 )L 2 (J^ J 2 ) (32) 


where 


L 0 (J*,J) s ^-± P 


6(J’,J) 


(33) 


L 2 (J’,J) = - ( J + _ !)( J + j +2) + - (J + !)( 2J + l )jL 6 (j« > j) + 3 ( J - _ 1 )J 5(j> j _ 2) 

4 2J + 3 2 (2J + 3)(2J - 1) 4 2J - 1 


(34) 
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Final Integrated Absorption Coefficient and Relation 
to Other Absorption Coefficients 


The final equation for the integrated absorption coefficient s^, 
stituting equation (32) into equation (9) 


is obtained by sub- 


s 


U' 


v 2 n 2 P(v v J v 0)P(v 2 ,3 2 ,0)W 

3hce o 


(35) 


where 


w . L 0 (i{,J 1 )i^(ij,jj)G 1 + l 2 (/j;,j ] )l 0 (.:j,;,j 2 )g 2 + l 0 (j;,j 1 )l 2 (j^,j 2 )g 3 

+ L 2 (J;,J 1 )L 2 (J^,J 2 )G 45 (36) 

G 1 ' ( a l«2 - A 2 ) 2l l < 37 > 

G 2 ■ Q 1 q2i 6 + i l Q 2 I 2 * a ^3 + *>1 a l<$7 - 2 <3lQ 2 a 2"8 ‘ 2 «1 a l‘3 2 a 2 ! 9 < 38 > 

G 3 ‘ «? Q 2 I 6 + 4 1 Q 2 I 3 + a ^2 - a l«2 I 8 + ^2 a 2 I 7 * a l«2 a 2 J 9 < 39 > 

G 45 " ( a l«2 - «1 a 2 ) 2<1 4 + 2I 5> (40) 

I. = E,(R)4»R 2 dR (1=1, 9) (41) 

Jo 


E i D 0000’ D 2000’ D 0020’ D 2020’ D 2-222» C 2000’ C 2000 D 2000’ C 2000 D 0020’ and 
D 2 qoo D 0020 for 1 = 1 t0 9 ’ respectively. 

In applying equations (35) to (41) there are several subtleties to note. There is no 
restriction on Vj or V 2 , so excited vibrational states may be included. There is no 
restriction on which v changes or how much it changes. In fact, one or both may 
change in a transition. Also, there is no restriction on whether v^ or v 2 increases or 
decreases provided that the wave number of the resulting transition is positive. Neither, 
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one or both J's may change, but the changes are limited to J' - J = ±2. If equation (35) 
is applied to a fundamental vibrational transition, the transitions with changing are 
considered distinct from those with Vg changing. 

The effects of mechanical anharmonicity and vibration- rotation interaction are con- 
tained in equations (35) to (41). Specifically, these effects influence the values of 
P( v l’Jl’0)’ Qi> Q2 > a i> and a 2‘ The model used in deriving equations 

(35) to (41) is called the vibrating- rotator model for brevity. 

The integrated absorption coefficient s of a transition with the wave number re- 
moved was defined by equation (3) but is density dependent. To remove this density de- 
pendence we define a new integrated absorption coefficient S for the £ — £' transition. 



2 

P 


where p is dimensionless Hg density given by 


n 

P H — 


(42) 


(43) 


and n Q is the Loschmidt number. If s^, is appropriately summed, equation (42) can 
be modified so S applies to a rotational branch or branches or to all transitions with a 
net change of +1 vibrational quanta. 

The spectral absorption coefficient with the wave number removed is defined by 


A 


v 



Then 


S = 



oo _ 

A v dp 


(44) 


(45) 


Assuming local thermodynamic equilibrium, the spectral absorption coefficient in- 
cluding wave number and stimulated emission is (ref. 21) 


A- = A~v [ 1 - e 

v v 


-hcu/kT 


( 46 ) 
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so the attenuation of a monochromatic light beam by a path length l in gas is 


a 2 . 

I - A v p 1 

— = e 


(47) 


Line Shapes 

The line shape for a transition depends on what mechanism induces the electric di- 
pole moment that makes the transition possible. For pressure-induced transitions there 
are two mechanisms: (1) the quadrupole moment of a molecule can induce a dipole 
moment in another molecule or (2) the overlap of the electron orbitals of two molecules 
can distort the orbitals and thereby cause a dipole moment. For fundamental vibrational 
transitions at room temperature and below, the principal mechanism has been found to 
depend on which rotational branch the transition belongs to (ref. 11). Only transitions 
in which the total angular momentum of not more than one molecule changes are impor- 
tant. Such transitions are classified into O, Q, or S branches according to whether J 
changes by -2, 0, or +2, respectively, in an absorption transition. 

At room temperature the Q branch is principally due to the overlap dipole moment 
whereas the O and S branches are principally due to the quadrupole- induced dipole 
moment. Poll (ref. 21) has stated that line profiles in vibrational and rotational spectra 
which are due to the same quadrupole-induced effect have the same shapes, and those due 
to overlap effect have about the same shapes. At higher temperatures it is impossible to 
separate the quadrupole and overlap effects, so we assume one shape for all Q-branch 
lines and another shape for all O- and S-branch lines. The shapes are based primarily 
on experimental results. For both shapes 


A~ = 


Sf~ 


v rxi 

/ f ~ 

J 0 v 


du 


(48) 


where f~ is the line shape function. 

Q-branch lines. - From Hunt and Welsh (ref. 11) 



w q exp 


• he 


V - u' 

o 

kT 


~ ,2 2 

(v - u Q ) + w ( 


Q 



(49) 
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(50) 


w, 


f~ - 

V 


Q 


& 


~ x 2 


(u Q =s P) 


+ w, 


Q 


where is the wave number of the line computed from the term values, and Wq is 
the half width given by 


/ T _ 1 

w„ = 290( — — ) cm 1 (51) 

Q \300°k) 

Watanabe and Welsh (ref. 22) noted that for H 2 - He mixtures the high wave number wing 
was exponential for (v - v )/w.-. > 0. 7. In addition, for theoretical reasons all moments 

O 

of a line should be finite (private communication from V. F. Sears, Atomic Energy of 
Canada, Ltd. , Chalk River, Ontario). This requirement is satisfied by an exponential 
tail but not by a dispersion tail. Hence we assume 


f~ = F exp 



(V Q + 0. 7wq < V < 00 ) 


(52) 


where F and U are constants chosen so equations (50) and (52) agree in value and 
slope for v - v Q + 0. 7 Wq. 

O- and S-branch lines . - From Mac Taggart and Hunt (ref. 23) the shape of pure ro- 
tational lines of pure hydrogen is 



Wg exp 




(v 


V 


+ w 


2 

S 


(v 


V 


w c 

f~ = ± 

v 2 


{v-v Q ) +w s 


(v Q ^ V < v Q + 1. 73wg) 


(53) 


(54) 



M 


, o)2 .40 


(v + 1. 73w g < V < P Q + 7w g ) 


(55) 
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where the line width is 


Wg = 133^ ^-Y^cm”* (5 

b \298 K/ 

and M was chosen so equations (54) to (55) have the same value at u = + 1. 73wg. 

Because pure rotational lines are principally quadrupole induced, equations (53) to (56) 
can also be used for vibrational O and S lines. For the same reason as before, we 
assume an exponential tail. 


/ v ~ *'V\ 

f~ = G exp (-H -] (u + 7w g <= v < «>) 

" V W S ) ° 


(57) 


where G and H are constants chosen so equations (55) and (57) agree in value and slope 
for v = v Q + 7wg. 


RESULTS AND DISCUSSION 
Dipole Moment Integrals 

Calculating the dipole moment integrals 1^ to Ig requires the expansion coeffi- 
cients of the dipole moment and its derivative and also the inter molecular potential. The 
major contributions to the dipole moment integrals come from R = 2. 5 to 8 bohr 
(1. 323x10”*® to 4. 233x10” 1** m) so reliable expansion coefficients and potential are re- 
quired for this range. 

Expansion coefficients were obtained by five methods, the method used depending on 
the value of R. Expansion coefficients for R = 2. 5 to 4 bohr (1. 323x10”*® to 
2. 1 17x10” m) were obtained by an orthogonalized valence-bond calculation including 
full configuration interaction (ref. 12). The results are given in table I. For R = 6 
to 8 bohr (3. 175x10”*® to 4. 233x10”*® m) they were obtained by the large-R approxima- 
tion (ref. 13). Both sets are shown in figures 3 and 4. For R = 4 to 6 bohr (2. 117x10”*® 

to 3. 175x10”*® m) the two sets were faired together by means of curves of the form: 

4 2 

Coefficient = c jR“ + Cg + c^R + c^R , where c^ to c^ are constants. This is also 
shown in figures 3 and 4. For R = 0 all the expansion coefficients must be zero. Hence 
for R = 0 to 2. 5 bohr (0 to 1. 323x10”*® m) the expansion coefficients were interpolated. 
For R > 8 bohr (4. 233x10”*® m) they were extrapolated. 

The inter molecular potential was obtained by one method for large R and another 
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for small R. For R > 4 bohr (2. 117x10”*® m) the Morse potential of Fisher (ref. 24) 
was used. Fisher fitted his potential to eight sets of experimental data (refs. 25 to 32). 
The viscosity data (refs. 25 to 31) were taken at temperatures between 200 and 2341 K. 
Second virial coefficient data (ref. 32) were also included by Fisher. For R < 4 bohr 
(2. 1 17x10“ m) an exponential was fit to the Morse potential at 4 bohr (2. 117x10”*® m) 
and the orthogonalized valence-bond potential (table I) at 2. 5 bohr (1. 323x10”*® m) (see 
fig. 5). 

The dipole moment integrals were integrated on a digital computer between R = 0 
and 25 bohr (0 and 13. 23x10"*® m) and are given in table II. 


Integrated Absorption Coefficients 

Equations (35) to (42) were used to calculate the integrated absorption coefficients S 
of transitions with a net change in vibrational quanta of +1. Vibration- rotation term 
values were obtained by adding rotational term values of Patch (ref. 33) to vibrational 
term values of Spindler (ref. 34). Only states with energies below the dissociation 
energy were included. Vibrating- rotator vibrational overlap integrals and matrix ele- 
ments were calculated on a digital computer using Spindler's potential energy to which 
2 2 2 

h J(J + 1)/8 tt m r r was added. Typical vibrational matrix elements are shown in fig- 
ures 6 and 7 and are compared with harmonic -oscillator rigid -rotator vibrational matrix 
elements used by all previous investigators. For large values of J the two models give 
very different values. Some vibrational overlap integrals are shown in figures 8 and 9 
and are compared with harmonic -oscillator rigid-rotator overlap integrals used by all 
previous investigators. Again there are some large differences. The vibrating- rotator 
matrix elements and overlap integrals agree well with James (ref. 35). It would be 
reasonable to expect that the large differences between the results from the two models 
would have large effects on S. 

Preliminary integrated absorption coefficients S for 298 and 5200 K obtained by 
evaluating equations (35) to (42) on a digital computer showed that certain transitions 
could be neglected. Double vibrational transitions where one molecule had a change of 
+2 vibrational quanta and the other had a change of - 1 vibrational quanta had negligible S 
values. Because of this and because all double vibrational transitions are forbidden for 
the harmonic- oscillator rigid-rotator model, all double vibrational transitions were 
henceforth neglected, and only single vibrational transitions where one molecule had a 
change of +1 vibrational quanta were included. In the latter, double rotational transi- 
tions had negligible S values and were therefore excluded in the final program. How- 
ever, the final program did include double transitions in which there was a vibrational 
transition in one molecule and a rotational transition in the other molecule. 
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Final integrated absorption coefficients summed over all values of Vj, Vg, Jj, and 
J g are given in figure 10. At all temperatures most of the contribution to the total inte- 
grated absorption coefficient is due to Q branches. 

A comparison of values of S for single vibrational transitions with a change in v 
of +1 for three models is given in figure 11. The vibrating- rotator model gave S as 
much as 30 percent higher than the simpler harmonic -oscillator rigid-rotator model. A 
hybrid model with vibrating-rotator term values and partition functions and with 
harmonic -oscillator rigid-rotator vibrational overlap integrals and matrix elements is 
also included. 


Spectral Absorption Coefficients 


The spectral absorption coefficient with wave number removed was calculated from 
S by equations (48) to (57) and is given in figure 12 for the vibrating-rotator model for 
single vibrational transitions with a change in v of +1. The major peaks are due to Q 
branches. The small bump in the 600 K curve at about 4700 reciprocal centimeters is 
due to S branches. The small bump in the 7000 K curve at about 500 reciprocal centi- 
meters is due to O branches. The major peak shifts to smaller wave numbers as tem- 
perature is increased because of the smaller wave numbers of Q-lines with excited ini- 
tial states, the increased importance of O-branches at high temperature, and the de- 
creased importance of S branches at high temperature. 

Once the absorption coefficient has been calculated accurately for a number of tem- 
peratures and wave numbers, it can be calculated much more rapidly for any tempera- 
ture or wave number in the same range by using an approximate fit to the accurate re- 
sults. The absorption coefficient A~ can be approximated by a modified Q-line shape. 


rwj exp £j (v - v c ) 

Av ~~ .2 2 

(v - v c ) + Wj 


& < V 


A u a 




(v - £ c ) 2 +W l 


( V c~ p - p c + l w 2) 


A~ ~ b exp 


- (* - g c > 

g 


( p c + l w 2 - v ) 


(58) 


(59) 


(60) 


19 


Illlllllllllllllllllllllllllllllllll 


where the wave number v of the line center is given by 

log 10 (4172 - v c ) = -5. 1972 + 2. 1 log 1() T (61) 

and 

y = 6. 0273x10" 10 + 2. 2905x10" 13 T + 4. 0848x10" 17 T 2 (62) 

w x = 363. 96 + 1. 3530T - 3. 5807X10" 4 T 2 + 3. 3618X10" 8 T 3 (63) 

log 10 j = 161. 45T" 1 - 2. 6996 - 1. 9537X10' 4 T (64) 

w 2 = - 108626T" 1 + 697. 59 + 0. 14353T (65) 

log 1() b = 28. 765T" 1 - 9. 0461 + 1. 1552X10" 4 T (66) 

log 1Q g = 1. 4860 + 0. 44462 log 1Q T (67) 


The constants in equations (61) to (67) were determined by fitting the absorption coeffi- 
cients in figure 12 and similar data for wave numbers from 8000 to 40 000 reciprocal 
centimeters and are valid for temperatures between 600 and 7000 K and wave numbers 
between 100 and 40 000 reciprocal centimeters. In equations (58) to (67) temperature is 
in degrees Kelvin, and wave number is in reciprocal centimeters. On a high-speed digi- 
tal computer the use of equations (58) to (67) instead of equations (33) to (40), (42), and 
(48) to (57) results in two orders of magnitude less running time, one-fifth the storage 
locations, and accuracy at any wave number within ±10 percent of the maximum value 
A~ for the given temperature. 

From equations (46) and (47) the attenuation of a monochromatic light beam by a path 
length l in Hg gas in terms of A~ is 


— = exp 
X o 


- e- hc ^ kT )p 2 l 


(68) 


where A~ and p are dimensionless. 

For quick reference A~ is more convenient that A~ and was obtained from A~ by 
equation (46). The results for the vibrating- rotator model are shown in figure 13. 
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The calculation of A~ and A~ is greatly complicated by using the vibrating- 
rotator model rather than the harmonic -oscillator rigid-rotator model, so it is interest- 
ing to see how this affects the answer. This is shown in figure 14 for 4000 K. A hybrid 
model is also included in figure 14. The principal effect of using the vibrating- rotator 
model rather than the harmonic- oscillator rigid-rotator model is to shift the absorption 
to smaller wave numbers. 


Comparison with Other Investigators 

Hare and Welsh (ref. 10) measured the binary and ternary integrated absorption 
coefficients for the Hg fundamental at 298 K. After allowing for their different definition 
of the integrated absorption coefficient, we find the vibrating- rotator model gives 87 per- 
cent of their value for the binary integrated absorption coefficient. The discrepancy is 
principally due to the neglect of x and y components of the dipole moment derivative 
in the vibrating- rotator model. 

Linsky (ref. 7) calculated A~ by extrapolating Van Kranendonk’s (ref. 6) theory. 

The vibrating- rotator model is compared with Linsky's model in figure 15 for 3000 K. 

The vibrating- rotator model (labelled "This report") has an integrated absorption coef- 
ficient 2. 1 times Linsky's value. However, this is not so much due to the vibrating- 
rotator model itself as to the inter molecular potential used with it (see fig. 5). Linsky 
followed Van Kranendonk (ref. 6) in using the unrealistic Lennard-J ones 6-12 potential 
also shown in figure 5. 

Linsky's pressure-induced translational, rotational, and vibrational absorption coef- 
ficients together with absorption coefficients for 12 other processes were used by Patch 
(ref. 36) to calculate the Planck and Rosseland mean opacities of hydrogen. The results 
are shown in figures 16 and 17. The discontinuities in the curves at 3000 K are because 
Linsky gave no vibrational absorption coefficients above this temperature. Patch's cal- 
culation (ref. 36) was repeated using equations (58) to (67) for pressure-induced vibra- 
tional absorption instead of Linsky's. The resulting opacities for pressures of 100 and 
1000 atmospheres (1.013x10^ and 1. 013x10** N/m'*) are labelled "This report” in figures 
16 and 17. Below 3000 K equations (58) to (67) resulted in Planck and Rosseland mean 
opacities as much as factors of 1. 9 and 303 higher, respectively, than in reference 36. 
Above 3000 K equations (58) to (67) resulted in Planck and Rosseland mean opacities as 
much as factors of 18. 9 and 6. 5 higher, respectively, than in reference 36. 


Effect on Gaseous-Core Nuclear Rockets 

The opacity of hydrogen (figs. 16 and 17) is still not enough below 6000 K to obtain 
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adequate heat transfer to hydrogen and to protect the walls of gaseous-core nuclear 
rockets. Some opaque additive to the hydrogen would be necessary. 


CONCLUSIONS 

The coefficient for pressure- induced vibrational absorption in Hg - Hg collisions 
was calculated for temperatures between 298 and 7000 K and wave numbers from 100 to 
40 000 reciprocal centimeters with many refinements. The conclusions were 

1. The integrated absorption coefficient at 298 K was 87 percent of the experimental 
value. 

2. At 3000 K the integrated absorption coefficient was 2. 1 times a previous theoreti- 
cal estimate, principally due to use of a more realistic inter molecular potential in this 
report. 

3. Although the vibrating- rotator model gave quite different vibrational matrix ele- 
ments and overlap integrals compared to the usual harmonic -oscillator rigid-rotator 
model, it did not cause more than a 30-percent change in the integrated absorption coef- 
ficient. 

4. The vibrating- rotator model shifted the absorption to smaller wave numbers than 
for the harmonic-oscillator rigid-rotator model. 

5. Below 3000 K, where a previous estimate of vibrational absorption was available, 
absorption coefficients from this report resulted in Planck and Rosseland mean opacities 
as much as factors of 1. 9 and 303, higher, respectively, than the best previous esti- 
mates for pressures from 100 to 1000 atmospheres (1. 013x10^ to 1. 013x10** N/m 2 ). 

6. Above 3000 K, where no previous estimate of vibrational absorption was avail- 
able, absorption coefficients from this report resulted in Planck and Rosseland mean 

opacities as much as factors of 18. 9 and 6. 5 higher, respectively, than the best pre- 

7 8 

vious estimates for pressures from 100 to 1000 atmospheres (1.013x10 to 1.013x10 
N/m 2 ). 

7. The opacity of hydrogen below 6000 K is still not enough to obtain adequate heat 
transfer to hydrogen and to protect the walls of gaseous-core nuclear rockets. Some 
opaque additive to the hydrogen would be necessary. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, October 20, 1970, 

122-28. 
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APPENDIX - SYMBOLS 


V k) 

a~ 

v 

B 

b 

^ 1 Z 2^2 

c 

c t (i=l, 2, 3, 4) 

° Z 1 ^ 1 Z 2^2 
E i (i=l, 9) 

e 

F 

G 

G.(i=l, 2,3,45) 

g 

H 

h 

I.(i=l,2,...,9) 

I 


spectral absorption coefficient including wave number and stimulated 
emission 

spectral absorption coefficient with wave number removed and ex- 
cluding stimulated emission 

function of molecule k involving integral of normalized associated 
Legendre functions 

linear absorption coefficient for wave number v excluding stimu- 
lated emission 

integrated Einstein coefficient of absorption 

spectral Einstein coefficient of absorption 

function of temperature in modified Q-line shape 

expansion coefficient for p° 

speed of light 

constants for faired curve 

expansion coefficient for p ^ 

product of two expansion coefficients 
charge of electron 
constant in Q-branch line shape 
line shape function 

constant in O- and S- branch line shape 
functions of Q, A, and I 

function of temperature in modified Q-line shape 

constant in O- and S-branch line shape 

Planck constant 

dipole moment integral 

transmitted light intensity 

incident light intensity 



2) 


J!(i=l, 2) 


3 

k 

L.(i=0, 2) 

l 


M 

M.(i=l,2,...,5) 
m.(i=l, 2) 


m!(i=l,2) 


N 

n 


P < v i’ J i> m i>l 
(i=l,2) J 

P 

Q k (k=l, 2) 

R 

r k (k=l,2) 

S 


s 

T 

t 

U 

V 


total angular momentum quantum number for initial state of mole- 
cule i 

total angular momentum quantum number for final state of mole- 
cule i 

function of temperature in modified Q-line shape 
Boltzmann constant 
functions of J' and J 
path length 

constant in O- and S-branch line shape 
functions of Q, A, C, and D 

quantum number for component of total angular momentum along z 
axis for initial state of molecule i 

quantum number for component of total angular momentum along z 
axis for final state of molecule i 

reduced mass of oscillator 

number of Hg molecules 

number density of Hg molecules 

Loschmidt number 

probability of molecule i being in state with quantum numbers having 
values v., rm 

number of pairs of molecules 
vibrational overlap integral for molecule k 
intermolecular distance 
internuclear distance of molecule k 

integrated absorption coefficient with wave number and density 
removed 

integrated absorption coefficient with wave number removed 

temperature 

time 

constant in Q-branch line shape 
volume 
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V|(i=l, 2) 
v!(i=l, 2) 
W 


W Q’ W S’ W 1’ W 2 

x,y,z 

x i 

Y Jm 


Z 3 

z i 


r 

A k (k=l, 2) 
6(J',J) 


0 Jm 

^(i=l, 2) 


My? M z 

Mi(i=l, 2 ) 





o’ 


V 


c 


p 


dr 

$ 


vibrational quantum number for initial state of molecule i 
vibrational quantum number for final state of molecule i 
transition-dependent factor in s 
line widths 

Cartesian coordinates fixed in space (fig. 1) 
position coordinates of electrons 
spherical harmonic 

th 

number of elemental charges on the j nucleus 

z coordinate of i in electron 
th 

z coordinate of j nucleus 

function of temperature in modified Q-line shape 

vibrational matrix element for molecule k 

Kronecker delta function 

electric permittivity of free space 

normalized associated Legendre function 

polar angle of molecule i 

components of electric dipole moment 

3ju. z /3r^ evaluated at r^ = r° and rg = r^ 

\i for configuration 4 
z 


components of electric dipole moment matrix elements 


photon wave number 

wave number of line center 

dimensionless density (see eq. (43)) 

radiation energy density per unit wave number 

element of volume in configuration space 

inter molecular potential averaged over orientations 



<Pj(i-l> 2) 

* 

Vk (k=1,2) 

Subscripts: 

e 

n 

t 

V 

P 

P' 

Superscripts: 

o 

* 


azimuthal angle of molecule i 
wave function 

vibrational wave function of molecule k 

electronic 

nuclear 

v r v 2 , Jp J 2 collectively 
VpVjpJpJ^ collectively 
v p v 2 ,Jp J 2 , mp m 2 collectively 
v l’ v 2’ J i J 2’ m i m 2 c °ll ect ively 

equilibrium internuclear distance of diatomic molecule 
complex conjugate 
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TABLE I. - DIPOLE MOMENT EXPANSION COEFFICIENTS AND INTERMOLECULAR 


Inter - 
molecular 
distance, 

R, 

m 


POTENTIAL FROM ORTHOGONAL IZ ED VALENCE 

Expansion coefficient 


u 2000’ 
C m 


BOND CALCULATIONS 


1. 323x10" 10 
1. 588 
1.852 
2.117 


1 . 758 X 10" 30 
9 . 510 X 10 " 31 

4 . 947 
2.498 


D 0000’ 

C 


9. 806X10" 20 

5. 609 

2. 986 
1. 516 


D 2000’ 

C 


6 . 351 X 10*" 20 

3 . 709 
1 . 996 
1.024 


D 0020’ 

C 


- 0 . 169 X 10" 20 
169 
-. 172 
-. 130 


^ 2020 ’ 

C 

-0. 121X10" 20 
.009 
.004 
~. 006 


D 

2 - 222 ’ 

C 

0 . 077 X 10" 20 
.045 
. 027 
.016 


Inter - 
molecular 
potential, 

*, 

j 

6. 210x10" 19 
2 . 645 
1 . 127 

4 . 732 X 10" 20 


TABLE II. - DIPOLE MOMENT INTEGRALS 


! Temper- 
ature, 
T, 

h> 

I 2> 

I 3’ 

Dipo 

h’ 

le moment ini 
I 5’ 

tegral 

I 6’ 

i 7 , 

X 8’ 


K 

C 2 m 3 

C 2 m 3 

C 2 m 3 

C 2 m 3 

C 2 m 3 

C 2 m 5 

C 2 m 4 

c 2 m 4 

C 2 m 3 

600 

9. 980X10 -70 

4. 845x10" 70 

4.781X10' 71 

2. 339X10' 73 

5. 534x10 74 

4. 368X10' 91 

1. 415xl0" 80 

-4. 418X10' 81 

-1. 349x10 

1000 

1. 543x10" 69 

7.325 

5.363 

2. 550 

1.085X10" 73 

5.780 

2. 007 

-5. 313 

-1. 727 

2000 

2.991 

1. 389X10" 69 

6. 519 

2.889 

2.498 

9. 653 

3.603 

-7. 370 

-2. 580 

3000 

4. 515 

2.076 

7. 397 

3. 123 

3.891 

1. 380X10" 90 

5. 291 

-9. 186 

-3. 323 

4000 

6.074 

2. 774 

8.090 

3. 316 

5. 230 

1.809 

7.022 

-1. 080X10" 80 

-3. 978 

5000 

7. 648 

3.475 

8.656 

3.494 

6. 511 

2. 247 

8.773 

-1. 225 

-4. 564 

6000 

9. 227 

4. 176 

9. 136 

3. 674 

7. 747 

2. 688 

1. 053X10" 79 

-1. 357 

-5. 096 

7000 

1. 079X10" 68 

4.868 

9. 542 

3.876 

8.938 

3. 131 

1. 228 

j-1. 478 

-5. 583 
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Molecule 1 


Molecule 2 



Figure 1. - Cartesian and polar coordinates for collisions of two molecules. Protons are 
located at a, b, c, and d. The z axis passes through the midpoints of the molecules. The 
angles 0j and 0 2 are polar angles. The angles (p^ and ^ are azimuthal angles 
measured in the x-y plane. 




Molecule 1 Molecule 2 

Figure 2. - Configurations of two H 2 molecules for which |i°(R) and ^(R) 
are required by the theory in this report. 
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Expansion coefficient, Dqqqq and D 2 ooO' 


o 


10 


-29, — 




10 


-23 



1.5 


1 _ J_ 


2.0 2.5 3.0 3.5 

Inter molecular distance, R, m 


J 

4.0x10 


-10 


Figure 3. - The three largest expansion coefficients of the dipole moment and its derivative. 
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.008 — Orthogonalized valence-bond 

calculation (table I) 

\ Large R approximation (ref. 13) 



Intermolecular distance, R, bohr 


1.5 2.0 2.5 3.0 3.5 4.0xl0~ 10 

Intermolecular distance, R, m 

Figure 4. - The three smallest and least important expansion coefficients of the 
derivative of the dipole moment. The oscillation of D 2 o 2 q results partly from 
the interpolation formula used between the four points from orthogonalized 
valence-bond calculations and may not be entirely valid. However, D^q makes 
a negligible contribution to the pressure-induced vibrational absorption coeffi- 
cient. 


co 

co 



Vibrating rotator U 1 - J as noted) 


2 SxlO’ 11 . 5r — Harmonic oscillator rigid rotator (all J 1 - J) 



Total angular momentum quantum number, J 


Figure 6. - Vibrational matrix elements for H 2 for v = 0 and v' - 1 for two models show 
the effect of vibration -rotation interaction. 


Vibrating rotator (J* - J as noted) 

Harmonic oscillator rigid rotator (all J' - J) •*’ ' ' 

yO 



Total angular momentum quantum number, J 


Figure 7. - Vibrational matrix elements for H 2 for v=0 and v' = 0 for two models show 
the effect of mechanical anharmonicity and vibration rotation interaction. 
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Vibrational overlap integral, Q Vibrational overlap integral, 





Integrated absorption coefficient, S, 1/cm 


GO 
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Term values and Vibrational overlap 

partition integrals and 

functions matrix elements 

Harmonic oscillator Harmonic oscillator 

rigid rotator rigid rotator 

Vibrating rotator Harmonic oscillator 

rigid rotator 

Vibrating rotator Vibrating rotator 




Temperature, T, K 

Figure 10. - Integrated absorption coefficient for various branches of all single 
vibrational transitions with a change in v of +1. The total for all branches 
is also given. The vibrating -rotator model was used. 



Temperature, T, K 


Figure 11. - Integrated absorption coefficient for all single vibrational transi- 
tions with a change in v of +1, using three models. The and 

curves coincide for temperatures below 3000 K. 



Absorption coefficient with wave number removed, 

Absorption coefficient, A~, 1/cm dimensionless 


I 


5X10" 9 



Figure 12. - Calculated pressure-induced vibrational absorption coefficient for all single 
vibrational transitions with a change in v of +1. The wave number has been removed 
from the absorption coefficient in the manner customary for pressure-induced processes. 
The effect of stimulated emission is not included. Temperature is given in thousands of 
degrees Kelvin. 


10x10 ' 6 



Figure 13. - Calculated pressure-induced vibrational absorption coefficient for all single 
vibrational transitions with a change in v of +1. The effect of stimulated emission is 
included. Temperature is given in thousands of degrees Kelvin. The attenuation of a 
monochromatic light beam may be calculated from I/I 0 = exp(-A~ p 2 l), where p is dimen- 
sionless H 2 density (see eq. 43). 
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Term values and Vibrational overlap 
partition integrals and 

functions matrix elements 


Harmonic oscillator Harmonic oscillator 
rigid rotator rigid rotator 



Figure 14. - Pressure-induced vibrational absorption coefficient calculated using 
three models. Temperature is 4000 K. All single vibrational transitions with 
a change in v of +1 were included. 
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I 



Absorption coefficient, A^j, 1/cm 


6xl0~ 6 



0 1 2 3 4 5 6 7 8xl0 3 


Wave number, v , cm' 1 

Figure 15. - Comparison of pressure-induced vibrational absorption coefficient in the funda- 
mental region according to two different theories. The temperature is 3000 K. 



Figure 16. - Effect of H? - H 2 pressure-induced vibra- 
tional absorption coefficients of this report on the 
Planck mean opacity of hydrogen. The discontinuities in 
the dashed curves at 3000 K are due to reference 7 giving 
no vibrational absorption coefficients above 3000 K. The 
discontinuities in all curves at 4000 K are due to refer- 
ence 7 giving no rotational absorption coefficients 
above 4000 K. 
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Pressure, 
atm (Him 2 ) 



Figure 17. - Effect of H 2 - H 2 pressure-induced vibrational 
absorption coefficients of this report on the Rosseland 
mean opacity of hydrogen. Discontinuities in dashed 
curves at 3000 and 4000 K result from reference 7 giving 
no vibrational and rotational absorption coefficients, 
respectively, above these temperatures. 
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NASA- Langley, 1971 24 E-5843 



